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ABSTRACT 

Aims. Investigations of stellar encounters in cluster environments have demonstrated their potential influence on the mass and angular 
momentum of protoplanetary discs around young stars. In this study it is investigated in how far the initial surface density in the disc 
surrounding a young star influences the outcome of an encounter. 

Methods. The numerical method applied here allows to determine the mass and angular momentum losses in an encounter for any 
initial disc-mass distribution. Based on a power-law ansatz for the surface density, oc r~ p , a parameter study of star-disc encounters 
with different initial disc-mass distributions has been performed using N-body simulations. 

Results. It is demonstrated that the shape of the disc-mass distribution has a significant impact on the quantity of the disc-mass and 
angular momentum losses in star-disc encounters. Most sensitive are the results where the outer parts of the disc are perturbed by 
high-mass stars. By contrast, disc-penetrating encounters lead more or less independently of the disc-mass distribution always to large 
losses. However, maximum losses are generally obtained for initially flat distributed disc material. Based on the parameter study a fit 
formula is derived, describing the relative mass and angular momentum loss dependent on the initial disc-mass distribution index p. 
Generally encounters lead to a steepening of the density profile of the disc. The resulting profiles can have a r~ 2 -dependence or even 
steeper independent of the initial distribution of the disc material. 

Conclusions. From observations the initial density distribution in discs remains unconstrained, so the here demonstrated strong 
dependence on the initial density distribution might require a revision of the effect of encounters in young stellar clusters. The steep 
surface density distributions induced by some encounters might be the prerequisite to form planetary systems similar to our own solar 
system. 

Key words. Methods: numerical - Protoplanetary discs - (Stars:) circumstellar matter 



1. Introduction 



There is increasing observational evidence that most, if not all, 
stars are initially surrounded by a circumstellar disc. For exam- 
ple, [LadireFl^ (|2000 ) found by examining the L-band excess 
of young stars in the Trapezium cluster a fraction of 80 - 85% 
of all stars to be surrounded by discs. With time these proto- 
planetary discs become depleted of gas and dust and eventually 
disappear (JHaisch et al. 2001 Hillenbrand|2002 Sicilia-Aguilar 



et al.|2006 [Hernandez et al |20071|Currie et al.|26081|Hernandez 
et al |2008[rMamajek|2009HMassi et al.|2010| >. Currently it is un- 
clear which of a variety of physical mechanisms dominates this 
evolutionary processes, among them there are internal processes 
like viscous torques (e.g. |Shu et al.||1987| l, turbul ent effects 
Klahr & Bod enheimer||2003| l and magnetic fields (Balbus & 



iawley|2002|i, while photoevaporation (|Scally & Clarke|2 001 ; 



Clarke et al.|2001||M atsuyam a et al.|2003||Jo hnstone et al.|2004| 
Alexander et al.|2005| |2006UErcolano et al.|2008| JDrake et~al . 
2009 ;|Gorti & Hollenbach|2009|l and tidal torques j[rleller|1993 



Clarke & Pringle 1993; Ostriker 1994; Heller 1995; Hall et al 
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are candidates for external disc destruction processes. 



The focus of the present numerical investigation is on the 
effect of gravitational star-disc interactions on the disc-mass dis- 
tribution and, therefore, the mass and angular momentum of the 
discs. A star-disc encounteiQ can cause matter to become un- 
bound, be captured by the perturbing star or pushed inwards and 
potentially be accreted by the central star. The extent to which 
this happens depends on the periastron distance, the mass ra- 
tio of the stars, the eccentricity and, moreover, the initial (pre- 
encounter) mass distribution in the disc. So far, limitations in 
observations of protoplanetary discs prevented to reveal consis- 
tently how the surface density of low-mass discs develops. This 
means that due to observational constraints there does not exist a 
unique predetermined initial state for the disc-mass distribution 
before an encounter. A wide variety of surface densities have 
been derived from fitting resolved millimeter continuum or line 



emission data with parametric disc structure models (e.g. Mundy 
et al. || 1996] |Lay et al.]|1997| ) or in co mbination with broa dband 
spectral energy distributions (S EP) fWilner et al.||2000| |Testi 
et al.|20011|Akeson et al.|20021|Kitamura et al.|2002[|Andrews 
& Williams 2007b ). While those studies have profoundly shaped 



1 The term 'star-disc encounter' means here encounters in which only 
one of the stars is surrounded by a disc, in contrast to disc-disc encoun- 
ters, which denominate encounters in which both stars are surrounded 
by discs. 
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our knowledge of disc structures, they are all fundamentally lim- 
ited by low angular resolution. 

The standard fitting method underlies the assumption that 
the surface density 2 has a simple power-law dependence of the 
form 



Table 1. Initial conditions for all simulations 



2(r) oc r 



(1) 



out to some cut-off radius (e.g. |Andrews & Wi lliams 2 007 a) . 
Estimates based on numerical power-law models fitted to obser- 
vational data lead to distribution indices p ranging roughly from 
p — to p = 2. Recent works even found anomalous results, e.g. 
|Isella et aT] ( |2009) l observed distribution indices of p < 0. 

Analytical approaches suggest different disc-mass distribu- 
tion indices, too. The most widely used model is that for a steady 
state viscous accretion disc with a surface density distribution 
index of p — 1 (e.g. Hartmann et al.| [l998). However, simulat- 
ing the evolution of protostellar discs that form self-consistently 
from the collapse of a molecular cloud core lead to a surface dis- 
tribut i on index of p = 1. 5 (Lin & Pringlej 1 990| |Hueso & Guillot 



2005 Vorobyov & Basu 2007J while treating magnetized disc 
material results in a flatter disc-mass distribution of p — 0.75 
( Shu et al.|2007| l. Taking this huge variety of distributions into 
account, one has to consider different initial disc-mass distribu- 
tions to evaluate their effect in star-disc encounters. 

Most previous numerical studies of star-disc encounters used 
only a single density distribution, mainly focusing on the case of 
a theoretically motivated r~ l disc-mass distribution (Hall et al. 



T996l|Hall|1997HPfalzner|20041|Olczak et al.|20061|Moeckel & 
Bally 2006; Pfalzner & Olczak 2007|. Star-disc encounters with 
different initial disc-mass distributions were only considered in 
a very limited way. Heller] ( fT995| ) performed numerical simula- 
tions of two different mass distributions (p — and p = 1.5) con- 
centrating on parabolic encounters with equal mass stars while 
|Hall| ( |1997| l investigated initial surface distributions of p = 
and p = 1 for close and penetrating encounters of an unity mass 
ratio. A study of a wide parameter range focusing on multiple 
initial disc-mass distributions is still lacking. 

Nevertheless, numerical studies of star-disc encounters only 
allow a sectional view on the processes since it is not possible to 
simulate each combination of encounter parameters. Earlier an- 
alytical studies by Ostriker ( 1994 ) did not suffer from this short- 
coming. In her study a first order approximation of the angular 
momentum loss dependent on the initial disc-mass distribution is 
given. However, the validity of her results is limited to large pe- 
riastron radii (for example r per i/rdi sc > 3 for MilM\ = 1), where 
the angular momentum loss is usually well below 10 %. Close or 
even penetrating encounters cannot be interpreted by this linear 
perturbation theory ( Ostriker|1994| |Pfalzner et al.|200 5b) mak- 
ing numerical studies indispensable in this regime. 

In this work the effects of star-disc encounters have been in- 
vestigated for a large parameter space considering most config- 
urations that can be expected in a typical young cluster. The in- 
vestigated mass distributions cover the entire range of the so far 
observed disc-mass distributions. 

Section 2 shortly describes the numerical methods and pa- 
rameter range used in this study while Section 3 presents the re- 
sults of the simulations including a fit formula for the mass and 
angular momentum loss depending on the initial disc-mass dis- 
tribution index p. A summary and discussion is given in Section 
4. 



variable 


value 


Outer disc radius, r disc 


100 AU 


Inner disc radius, r^asc 


10 AU 


Number of particles, yip^ 




Disc-mass, m^sc 


10- 4 M o (10- 3 M o ) 


Particle distribution index, b 


0,(1) 


Mass distribution index, p 


i 1 i 

2' L > 4 


Rel. perturber mass, M 2 /M t 


0.1-500 


Rel. periastron distance, r per i/rdi sc 


0.1-20 


Eccentricity, e 


1 


Simulation time, f cn d 


~ 3000 yrs 



2. Methods 

The encounter between a disc-surrounded star with a secondary 
star is modeled using a code that is based on the numerical 
method described in |Pfalzner| ( 2003) >. In our simulations only 
one star is initially surrounded by a disc. However, previous in- 
vestigations showed that star-disc encounters can be generalized 
to disc-disc encounters as long as there is no significant mass 
exchange between the discs (Pfal zner et al.|2005a| l. In case of a 
mass exchange the discs can be to some extent replenished, so 
that for very close encounters the mass loss determined in this 
study would be slightly overestimated. 

A summary of the disc properties can be found in Table 1. 
The disc is represented by 10.000 pseudo-particles, distributed 
according to a given particle distribution oc r~ b . Choosing this 
relatively low number of simulation particles is motivated by the 
aim to cover a large encounter parameter space. However, per- 
forming test simulations with 50.000 particles, shows that the 
lower resolution actually suffice for the here investigated prop- 
erties. 

The simulation particles move initially on Keplerian orbits 
around the central star. The disc extends from an inner gap of 
10 AU to 100 AU. The inner cut-off avoids additional complex 
calculations of direct star-disc interactions and saves computa- 
tion time. Any pseudo-particle that reaches a sphere of 1 AU 
around the central star is removed from the simulation and stated 
in a commonly used simplified approach as being accreted ( |Bate| 



et al.|2002| | Vorobyov & Basu|2005| |Pfalzner et aL|2008) . The 
vertical density distribution in the disc follows 



p(r, z) = po exp 



2H 2 (r) 



(2) 



where po is the unperturbed mid-plane particle density on the 
equatorial plane with po oc and H{r) the vertical half- 

thickness of the disc (see Pringle 1981 ). H{r) is chosen as H(r) = 
0.05 r according to a temperature profile of T — Tq ■ r _1 with 
T = 20K at the inner edge of the disc. 

The temporal development during the encounter is modeled 
using a fifth-order Runge-Kutta Cash-Karp integrator with an 
adaptive time step size control for the numerical calculations. 
Long-range interactions of the gravitational forces between the 
disc particles and the perturbing star are calculated using a hier- 
archical tree method (B arnes & Hut|1986"l ). 

The total simulation time of t en d « 3000 yrs for each en- 
counter corresponds to three orbital periods of the outermost par- 
ticles in the investigated standard disc with size rdi sc = 100 AU 
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around a 1 M Q star. This time span was found to be adequate for 
the calculation of all relevant quantities. 

The discs obtained at the end of the simulations are still de- 
veloping in the sense that taking into account viscous forces 
the eccentric orbits of perturbed disc pseudo-particles would re- 
circularize on time scales probably in excess of 10 5 yr. However, 
bound and unbound material already can be clearly distinguished 
shortly after the encounter. 

In this study we are interested in protoplanetary discs, which 
are usually of low mass compared to protostellar discs (Bate 
|201 l| l. Here low-mass discs mean m<ji SC <K O.lMi, where Mi 
is the mass of the disc surrounded star. The actual used disc- 
mass value is l(T 4 Mi. However, we tested as well the case 
widisc = 10~ 3 Mi and found no differences in the result when 
including self-gravity in the simulations. In addition, we per- 
formed test simulations including pressure and viscous forces 
using a SPH code (for a description of the code see Pf alzner| 
(2003 )), but only found negligible differences (< 3%) in the re- 
sults. Therefore, in most of our simulations self-gravitation and 
viscous forces were neglected in favor of higher computational 
performance ( |Pfalzner|2003l[2004l l. 

The actual values of mass and angular momentum loss in- 
duced by the encounters are obtained by averaging over several 
simulations with different seeds for the initial particle distribu- 
tion. The errors in the mass and angular momentum losses are 
determined as the maximum differences between the simulation 
results. They typically lie in the range of 2 - 3 %. 

Besides these obvious performance benefits the possibility 
to neglect the self-gravitation and pressure and viscous forces 
allows to further optimize computation time by treating the disc 
particles as pseudo-particles without a fixed mass during the 
simulation. 

The standard method to generate the initial disc-mass distri- 
bution in the simulations is by assigning each pseudo-particle the 
same mass ([Boffin et al.| T998l|Pfalzner|2004[|Olc"zak et al.|2006| 



Pfalzner & Olczak 2007 ). This approach requires to perform the 



whole suite of simulations afresh for each variation of the ini- 
tial disc-mass distribution. By contrast, here the different disc- 
mass distributions are realized by using a fixed pseudo-particle 
distribution and assigning their masses after the simulation pro- 
cess according to the desired density distribution in the disc. The 
implementation of such a flexible numerical scheme allows us- 
ing one suite of simulations for any initial disc-mass distribu- 
tion. Note, that the post-processing of the particle mass is only 
possible because of the restriction to low-mass discs where self- 
gravitation and viscous forces can be neglected. 

The particle masses are assigned to each pseudo -particle in 
the diagnostic step to calculate the encounter-induced losses and 
final mass distributions dependent on the initial distribution. For 
the example of a power-law disc-mass distribution with index p 
the mass of a particle, ;«;, depends on its radial position in the 
disc r\, the total disc-mass m^sc, the number of pseudo-particles 
n par t, their radial positions rj and the underlying power-law par- 
ticle distribution with index b. Since the number of pseudo- 
particles is limited, the discrete particle masses m,- are obtained 
by 



-p+b 



■ «disc • r. 



-p+b 



(3) 



The so established independence of mass and particle distribu- 
tion allows to improve the resolution of our simulations by plac- 
ing most of the pseudo-particles initially where the interaction 
between the stars and disc is strongest. 



In most of our simulations a constant particle distribution 
(b = 0) is used. This is done because the effects on mass and 
angular momentum are largest at the outskirts of the disc. Using 
a constant pseudo-particle density (b - 0) throughout the disc, 
ensures a higher resolution in the outer parts even for steep mass 
density profiles compared to the standard approach. By contrast, 
if one would be interested in processes that mainly concern the 
inner parts of the disc, like for example accretion, it would be 
preferential to use not a constant particle distribution but one 
that guarantees a high resolution close to the star. 

We have tested the method of a posteriori mapping of the 
particle masses against representative test calculations for dif- 
ferent initial pseudo-particle distributions and against the stan- 
dard method. The differences in the results where in each case 
negligible and thus justify the generalization of our results for a 
constant particle distribution (b = 0). 

Apart from the density distribution the outcome of an en- 
counter depends on the relative periastron distance r per i/rdi sc and 
the mass ratio between the involved stars M2/M1, where M2 
is the perturber mass. The here investigated parameter space - 
for a summary see Table 1 - is chosen in such a way that it 
spans the entire range of encounters likely to occur in a typi- 
cal young cluster, like the Orion Nebula Cluster. Recent inves- 
tigations show that the ONC with its high central density might 
be typical for young clusters with a mass > 1000M G ( |Pfalzner| 
2009 ). The lower limit of the perturber mass ratio was chosen to 
be MijM\ = 0.1 as for smaller mass ratios the influence of the 
perturber becomes insignificant. The upper limit, M2IM1 - 500, 
is determined by the maximum possible mass ratio in the ONC, 
which is given by the hydrogen burning limit as lowest mass 
(0.08 Mq) and 40 M , the approximate mass of the #'C Ori sys- 
tem, which contains the most massive star in the ONC. The 
limit where perturbations at large distances become negligible 
depends on the perturber mass. The inner edge of the disc marks 
the lower value of r pel -i/rdi sc =0.1. 

3. Results 

Although in principle any mass distribution of circumstellar 
discs can be studied, exemplary four disc-mass distributions are 
investigated here in more detail. To cover the entire range of 
numerically and observationally determined disc-mass distribu- 
tions a constant mass distribution (p = 0), representing the lower 
boundary of expected distributions, and a p = 7/4 distribu- 
tion, providing an upper limit, are considered. Additionally, a 
p — 1 distribution is investigated for comparison to previous 
star-disc encounter results and a p = 1/2 distribution, which is 
in the range of analytical results for low-mass discs (Shak ura"&| 
|Sunyaev| 1 973 [ |Pring le 1981 ) and similar to results found in re- 
cent investigations considering magnetized material (see Sec.[T]i. 

The present investigation focuses on coplanar, prograde en- 
counters on parabolic orbits (<? = 1) (see example in Fig.^) . 
Previous studies have shown that even for clusters as dense as the 
inner ONC region, most encounters in star clusters are expected 
to be close to parabolic ( |Larson|[T990l |Ostriker|[T99"4"l |Qlczak| 
et al. 2010 ). Such parabolic encounters provide the strongest im- 
pact of the perturber on the disc since for higher eccentricities 
the perturber only interacts shortly with the star-disc system and 
is, therefore, unable to influence the disc significantly. The limi- 
tation to a certain orientation is more restricting, as inclined and 
retrograde encounters can lead to lower losses i n mass (|Heller 



T993l[Hall et al.|[l996l |Clarke & Pringle 1993; |Pfalzner et al. 



2005bjt. However, Pfalzner et al. ( 2005b ) showed that up to an 
inclination of 45° the mass loss of the disc is only slightly re- 
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duced in comparison to a coplanar orbit. This means that our 
results have to be regarded as an upper limit for encounters at a 
different inclination. 



3.1. Surface density distribution 

To have a better understanding of the re-distribution of the disc 
material during a star-disc encounter in dependence of the ini- 
tial disc-mass distribution we first focus on the evolution of 
the surface density profiles. The evolution of the surface den- 
sity in star-disc encounters in dependence of the periastron dis- 
tance for a constant initial disc-mass distribution (p = 0) and 
one with a steep initial distribution (p = 7/4) can be seen in 
Fig. [T] It shows the mass distributions before (solid line) and 
after a penetrating (r per i/rdi sc < 1, dashed-dotted line), graz- 
ing (r per j/rdi sc = 1, double -dotted line) and distant encounter 
('"peri /'"disc = 3, dashed line). Here, the relative perturber mass 
was chosen to be M2/M1 = 1. 
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Fig. 1. Initial (solid line) and final surface densities in case of 
a) initially constant distributed disc material of p = and b) 
a steep distribution of p = 7/4. In both plots non-penetrating 
OperiMisc = 3, dashed line), grazing (r per j/r disc = 1, double- 
dotted line) and penetrating encounters (r per i/rdi sc =0.1, dashed- 
dotted line) are plotted for a perturbing star of equal mass 
(M 2 /M, = 1). 



As known from previous investigations ( |Hall 1997 Larwood 
|1997| l an encounter reduces the density in the outer parts of the 



disc by transporting some part of the outer disc material in- 
wards and some part migrating outwards. The latter might be- 
come captured by the perturber and unbound if the encounter is 
strong enough. These effects become more pronounced for small 
'"pen/^disc and/or large M2/M1. However, the actual amounts de- 
pend on the initial disc-mass distribution. In an equal-mass dis- 
tant encounter (r per i/rdi sc = 3) most of the perturbed disc ma- 
terial is pushed inside the disc towards lower disc radii. In this 
case 12.5 % of the total disc-mass migrates inwards for an initial 
constant disc-mass distribution (p = 0), whereas it is only 6 % 
for an initial p = 7/4. 

The effect of different initial disc-mass distributions becomes 
even more obvious for closer encounters like grazing (double- 
dotted lines in Fig. [T]l or penetrating fly-bys (dashed-dotted lines 
in Fig. [TJ. Here, the migration process and also the differences 
between the initial distributions are dominated by material mov- 
ing outwards. In case of a grazing encounter and an initially con- 
stant disc-mass distribution ~ 60% of the disc-mass becomes 
unbound. However, in addition ~ 10% of the total disc-mass 
can be found outside the initial disc radius of 100 AU still bound 
to the central star. By contrast, for the steep p = 7/4-mass dis- 
tributions it only becomes ~ 30 % of the disc-mass unbound but 
again ~ 10% is bound but situated outside 100 AU. We con- 
clude that generally the outer disc material is separated from the 
disc for grazing encounters while material initially located inside 
the disc (r/rdi SC < 0.7) will be re-distributed but remains bound 
to the central star. Hence, prominent differences in encounter- 
induced disc losses for the investigated disc-mass distributions 
are expected for strong perturbations of the outer disc parts. In 
penetrating encounters part of the disc material is pushed fur- 
ther inside the disc resulting in an increased surface density in 
the inner disc regions. In extreme cases the disc loss can be in- 
creased to > 90 % so that the final disc structure can no longer 
be regarded as a disc. 

In nearly all cases a steepening of the density profile is the 
general effect of an encounter. The degree of steepening depends 
on the mass ratio, periastron distance and the initial disc-mass 
distribution. Fig.[3]shows that even initially flat distributed disc- 
material (p = 0) can be redistributed into a surface density pro- 
file steeper than p as 2. 

How is the vertical structure of the disc influenced by an en- 
counter? As only coplanar encounters are investigated here, the 
answer is: surprisingly little. Fig. [2j? shows the vertical particle 
profile of a strongly perturbed disc after a penetrating encounter 
of M2/M1 = 1. The main effect of the encounter is a decrease in 
the number of particles in the outer disc regions. In these parts 
of the disc the resolution is somewhat diminished, but due to 
our constant particle distribution (b = 0) still significantly larger 
than in previous investigations using r particle distributions in 
their simulations. 

In a few special cases where material is being pushed moder- 
ately inwards even partially negative distribution indices as low 
as p = — 1 can be the end product of an encounter (Fig. [4]). The 
effect is most prominent for the inner part of perturbed discs 
with initially constant distributed disc material. Nonetheless, the 
outer parts of the disc always represent a distribution with posi- 
tive index, p > 0. 

This could perhaps explain the negative indices observed for 
some discs. Observations are still limited to measuring only part 
of the radial extension of the entire disc. If this happens to be 
the range were negative indices prevail, one would wrongly ex- 
trapolate negative indices for the whole disc, whereas in fact the 
overall index of the disc would still be positive. 
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a) 



b) 




Fig. 2. The dashed line in a) shows the boundaries of the initial 
disc while the solid line indicates the trajectory of a grazing per- 
turber (r per i/rdi sc - 1) of equal mass (MijM\ = 1). Material that 
resides within the disc after the perturbation is marked as black 
squares while material that is in the end bound to the perturbing 
star, unbound or accreted is shown as gray dots. Note that the 
simulations were performed in three dimensions as can be seen 
in b). 



1- 




0.01 ■ 



0.001 



^ Iodise 

Fig. 3. Initial (p = 0, solid line) and final surface density after 
an encounter of M 2 IM\ = 5.0 and r pel -i/r^ sc = 0.7 (dotted line). 
The dashed straight line represents a slope of p = 5/2. 




r/r d isc 

Fig. 4. Initial (p = 0, solid line) and final surface density after 
an encounter of M 2 jM\ = 1.0 and r per i/rdi sc =2.0 (dotted line). 
The dashed straight line represents a slope of p = — 1. 



3.2. Relative disc-mass and angular momentum loss 
3.2.1. Relative mass loss 

Besides changing the shape of the disc-mass distributions, star- 
disc encounters can remove material from the disc. Fig.[5]shows 
the relative disc-mass loss Am re i = (m t =o - m tcM )/m, = o for the 
four initial disc-mass distribution indices p = 0, p = 1/2, p = 1 
and p = 7/4. Two representative cases of orbital parameters 
are singled out: Fig. [5J? shows the dependence on the perias- 
tron distance for an equal-mass encounter, whereas Fig. [5|? de- 
picts the dependence on the mass ratio for grazing encounters 

(>"peri/>"di S c = 1)- 

Qualitatively, the dependency of the relative disc-mass loss 
on the periastron distance follows previous results ( Pf alzner etal.| 
2005b{ [Olczak e t al. 2006 ). However, the absolute values change 
considerably for the different initial disc-mass distributions. The 
effect is largest for nearly grazing encounters where mainly the 
outer disc regions are affected. As to be expected from the re- 
spective fraction of material in the outer regions, maximum mass 
losses are obtained for initially constant disc-mass distributions 
while losses for the r _7 ^ 4 -distribution are lowest. The largest dif- 
ference in mass loss between the investigated disc-mass distri- 
butions in Fig.^ occurs for r per i/rdi sc = 0.9, whereas the f 1 ^- 
distribution has only a mass loss of 33 %, the constant mass dis- 
tribution has a mass loss of 64 %. 

The situation is somewhat different for very close penetrat- 
ing encounters (r per ;/?di sc < 0.3) where the discs are so strongly 
perturbed that the resulting structure can hardly be described as 
a disc. In this case the disc-mass loss seems relatively indepen- 
dent of the initial mass distribution (see Fig. [3}j). At the other 
end of the parameter space - i.e. at large relative periastron dis- 
tances - the mass loss becomes too small (Am re i < 10%) to infer 
any dependence on the initial distribution. 

Similarly, the dependence on the initial mass distribution 
is less pronounced for M 2 /M l < 0.3 and M 2 jM\ > 90 (see 
Fig. [5j? for a grazing encounter). So generally weak perturba- 
tions - whether distant or of low mass ratio M 2 IM\ - are not capa- 
ble to significantly influence the discs, while in the case of strong 
perturbations nearly the entire disc material is removed indepen- 
dent of the investigated disc-mass distributions. In both cases the 
mass loss does not depend strongly on the disc-mass distribution. 
By contrast, encounters of intermediate strength are most sensi- 
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Fig. 5. The relative disc-mass loss of a p = (solid line), p - 1/2 
(dashed-dotted line), p — 1 (dotted line) and p = 7/4 (double- 
dotted line) disc-mass distribution including all particles bound 
tighter to the central star than to the perturber and excluding un- 
bound and accreted particles. The data is plotted for a) differ- 
ent periastron distances and an equal-mass perturber and for b) 
different perturber mass ratios and r per i/rdi S c = 1. The vertical 
dashed line indicates the initial outer disc radius. 



tive to the disc-mass distribution. For the case shown in Fig. |5j7 
we find maximum differences of up to 35 % for M2/A/1 = 3. 

In case of high-mass ratios M2/M1 > 20 and certain non- 
penetrating periastron distances, which restrict the gravitational 
star-disc interactions to the outer disc parts, differences in mass 
loss of even up to 40 % can be inferred for the different initial 
disc-mass distributions (see online appendix A). 

3.2.2. Relative angular momentum loss 

The different initial disc-mass distributions do not only influ- 
ence the disc-mass loss but also the angular momentum loss. 
Fig. |6]3 shows the relative angular momentum loss A7 re i = 
(Jt=o - Jt end )/Jt=o as a function of the encounter distance for 
equal-mass encounters and the four different disc-mass distribu- 
tions used in this work. As expected from previous results (e.g. 
|Pfalzner & Olc zak 2007 ) the general trends in relative angular 
momentum loss are quite similar to that of the mass loss (com- 
pare Fig. [5] and Fig. [6]) with angular momentum losses being 
slightly higher than the disc-mass losses. 



While the mass losses for encounters of intermediate 
strength with r per j/r d i sc = 0.9 and M 2 /Mi = 1 are 64% for the 
constant mass distribution compared to 33 % for the 7/4-mass 
distribution as shown before, the corresponding angular momen- 
tum losses are 75 % and 50 % , respectively. Material migrating 
inwards or becoming unbound due to an encounter leads to a rel- 
ative angular momentum loss, while part of the disc-mass which 
is pushed outwards the initial disc-radius but remains bound to 
the central star increases the total angular momentum of the disc. 
In total, the dependence on the initial disc-mass distributions is 
less pronounced for the relative angular momentum loss than for 
the disc-mass loss. 
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Fig. 6. The relative angular momentum loss of a p — (solid 
line), p — 1/2 (dashed-dotted line), p — 1 (dotted line) and 
p - 7/4 (double-dotted line) disc-mass distribution including 
all particles bound tighter to the central star than to the per- 
turber and excluding unbound and accreted particles. The data 
is plotted for a) different periastron distances and an equal- 
mass perturber and for b) different perturber mass ratios and 
'"pen/^disc = 1. The vertical dashed line indicates the initial outer 
disc radius. 



However, although the differences are lower in maximum 
the influence of different initial disc-mass distributions on the 
angular momentum loss covers a significantly large parameter 
range (see also online appendix B). 
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Since mass and angular momentum losses are generally in- 
fluenced by perturbations of the outer disc parts an initially flat 
particle distribution with high resolution of the outer disc regions 
results in an increased accuracy of the disc losses. Nevertheless, 
disc losses obtained with initially steep particle distributions 
(p = 1) found in previous studies (Olczak et al. 2006 Pfalzner & 
|Olczak|2007[ ) are completely reproduced within the error range 
in the present study. 

3.2.3. Adapting a fit formula dependent on the initial 
disc-mass distribution 

The numerical results for the mass and angular momentum loss 
in this study cover a wide parameter range but, however, present 
only a discrete classification of the relative losses for the differ- 
ent initial disc-mass distributions. Analytical approaches are a 
possible option to avoid this disadvantage but are only valid for 
a very limited parameter space of distant encounters (Ostriker 
D'Onghia et al.|[2010)> . To obtain a general estimate 



4. Discussion and conclusion 



1994 



of the effect of arbitrary initial disc-mass distributions on the 
encounter-induced mass and angular momentum loss of proto- 
planetary discs we present a fit formula that is valid for any ini- 
tial disc-mass distribution, given that it can be expressed by a 
power law of the form E(r) oc r ~ p . For this purpose we extended 
the fit function for the relative mass loss by |Qlczak et al.| p006 ) 
(see Eq.4 therein) and for the relative angular momentum loss by 
Pfalzner & Olczak (2007) (see Eq.l therein), which are valid for 
a r~ l -distribution, towards arbitrary surface density distribution 
indices < p < 7/4: 



Am r ,i = 




A/ re , = 1.02 



x exp 



M 2 l l -p expQ?) \ 
M m + \ 5 100 / 



0.5-r o 



(5) 



where r p = r peri lr4 isc is the relative periastron distance. 

For most of the parameter space the adopted functions fit the 
data well inside the error range and extend Ostriker's analytical 
function for the angular momentum loss considerably. Larger de- 
viations of the fit functions from the simulated losses occur only 
for high encounter mass ratios of Mi/M\ > 20, which was also 
the case for the established fit functions by Pfalzner & Olczak 
( 2007jl and Olczak et al. ( 2006 1. In case of disc-penetrating orbits 
and low-mass ratios the disc losses can be affected by Lindblad 
and Corotation resonances that are located in the inner disc re- 
gions and cause moderate deviations of the fit function from the 
expected disc losses. 

The extended fit functions provide a significant improvement 
of previous analytical and numerical results. They cover a huge 
parameter space of the most reasonable initial disc-mass distri- 
butions and the most relevant orbital parameter ranges expected 
for interactions in star clusters of any age. 



Most stars are believed to form in clusters (Lad a & Lada|2003) 
and probably undergo at least one encounter closer than 1000 AU 
during the lifetime of their disc (~ 10 6 yrs) ( Scally & Clarke| 
2001 ). So encounters are likely to have a significant effect on the 
disc structure. The scope of this study was to determine which 
role the initial mass distribution in the disc plays in this context. 
A full parameter study of star-disc encounters, as expected to 
occur in young dense star clusters like the Orion Nebula Cluster, 
has been performed with different disc-mass distributions of the 
form 2(r) oc r~ p , p e [0, 1/2, 1,7/4], that cover the whole range 
of observed mass distributions in discs. The main results can be 
summarized as follows: 

1. The relative disc-mass loss among the different initial den- 
sity distributions differs by up to 40 % for the same type of 
encounter. The largest differences are associated with strong 
perturbations of the outer disc edge, i.e. a grazing encounter 
in case of equal stellar masses. 

2. Although higher amounts of relative angular momentum 
than disc-mass are lost due to a star-disc encounter, the de- 
pendence on the initial disc-mass distribution is less pro- 
nounced. Nevertheless, half of the parameter range where 
angular momentum loss occurs shows differences among the 
investigated distributions of more than 15 %. 

3. The disc-mass and angular momentum losses due to a 
parabolic encounter can be fitted by a function depending 
on the perturber mass ratio, the relative periastron distance 
and the index of the initial disc-mass distribution. 

4. A steepening of the surface density slope is a general effect 
of an encounter for any initial disc-mass distribution and can 
result in distribution indices of p > 2 even in case of initially 
flat distributions. 

In short, the intuitive result that the flatter the mass distri- 
bution the stronger the change of the disc-mass and angular mo- 
mentum due to encounters has been quantified in this study. This 
infers that the importance of encounters in young stellar clusters 
and their potential to trigger encounter-induced losses might be 
reconsidered depending on the dominant initial disc-mass distri- 
bution. 

Another consequence of such an encounter is a change of the 
surface density in the disc on short time scales . In this context, 
the potentially most significant result of this study is that close 
encounters can provoke density profiles steeper than 2(r) oc r~ 2 
independent of the initial disc-mass distributions. As such den- 
sity profiles are claimed to be the prerequisite to form a planetary 
system similar to the solar system, encounters in the early history 
of the solar system could have provided these conditions (Adams 
2010). Desch (2007) analytically inferred an initial surface den- 
sity profile of the solar nebula protoplanetary disc of roughly 
2(r) oc r~ 2 - 2 to form the present planets of the solar system. As a 
mechanism for the formation of such a steep disc-mass density 
profile he considered photoevaporation by an external massive 
star (see also |Mitchell & Stewart|20T0l ). |Kuchner| 42004] ) found 
similar results for extrasolar planetary systems suggesting large 
disc-mass distribution indices of p — 2.0 + 0.5. 

Our results emphasize that these steep density profiles do 
not have to exist ab initio or be formed by photoevaporation pro- 
cesses but that even discs of initially constant distributed mate- 
rial can fulfill the requirements for the formation of a solar sys- 
tem type planetary system in the inner disc regions after a close 
encounter. Additional evidence for such an encounter in case of 
the early solar system was given by distant solar system objects 
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on highly eccentric orbits, like the transneptunian object Sedna, 
and also the sharp outer edge of the solar system at ~ 50 AU 
from the Sun Jlda et al.||1999[|Allen et al.||2001| |Morbidelli &| 
|Levison|2004}|Kenyon & Bromley |2004| l. 

On the other extreme, encounters can not only lead to profile 
steepening but as well profile flattening. Even unexpected sur- 
face distribution profiles of p < as observed by |Isella et aT] 
(2009| can be explained by the influence of a perturbing star. 
Non-penetrating encounters of a star-disc system with initially 
flat distributed disc material easily lead to such density profiles 
of the inner disc regions. However, photoevaporation as another 
significant external impact on protoplanetary discs would have 
to be considered as an additional driver in these considerations. 

Consequently, if all young stars would start out with the same 
disc density structure, the influence of the cluster environment 
by means of encounters in the very early and dense phases of 
cluster evolution could account for the observed multitude of 
disc-mass density profiles. 
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